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ABSTRACT 

We present HST/WFPC2 images of the symbiotic star V1016 Cyg which, for the first time, 
directly explore the inner regions of the nebula to AU scales at optical wavelengths. They also 
suggest that the [O in] AA4959,5007 emission takes place in a bipolar outflow. We use these 
images to determine the position of the two stars and hence a projected binary separation of 
84 ± 2 AU (assuming a distance of 2 kpc) and a position angle of 143.5 ± 10°. Furthermore, 
we combine our images with VLA radio imaging which enables diagnostic tests to be applied 
and properties of the circumstellar nebula to be determined. Finally we have obtained archive 
HST/STIS spectra of V1016 Cyg with which we were able to spatially resolve the source at 
various positions in the nebula. This enabled discovery of the ultraviolet counterpart to the ~25 
arcsec extended emission found by Bang et al. (1992). 

Subject headings: binaries: symbiotic — stars: individual (V1016 Cyg) — circumstellar matter — radio 
continuum: stars 

1. Introduction estimates in the literature, ranging from 2.1 to 10 

kpc (Watson et al. 2000, and references therein). 
V1016 Cyg is a symbiotic star which falls into The binary period of the system has been sug _ 

the D(usty)-type sub-category, containing a Mira gested to be 80±25 years from spectropolarimetry 

variable as its late-type cool component and a (Schild & Schmid 1996) or more recently as ~ 15 

white dwarf as the hot component (e.g. Taranova years from optical photometry (Parimucha et al. 
& Yudin 1983). There are a number of distance 
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2000). Infrared photometry also shows a ~450 day 
periodicity which has been attributed to the pul- 
sations of the Mira (Harvey 1974). 

The system is well known for its nova-like op- 
tical outburst in 1965 and has since been dubbed 
a symbiotic nova (see e.g. Parimucha et al. 2000 
and references therein). The outburst is thought 
to be the result of a thermonuclear runaway on 
the surface of the white dwarf as it accretes mate- 
rial emitted by the Mira (Mikolajewska & Kenyon 
1992). As with classical novae, symbiotic novae 
may be accompanied by a radio outburst - while 
radio monitoring did not take place in 1965, V1016 
Cyg is certainly one of the brightest radio sources 
among the symbiotic systems (e.g. Belczyhski et 
al. 2000). 

The source has been imaged previously in the 
radio and emits a spectrum typical of thermal 
bremsstrahlung. While originally thought to be 
produced in an expanding shell as a result of the 
1965 outburst, Watson et al. (2000) used MER- 
LIN observations to show that it is more likely that 
the peaks in the radio maps are the result of col- 
liding winds between the two stellar components. 
Despite a number of similarities between V1016 
Cyg and another D-type symbiotic, HM Sagittae, 
V1016 Cyg does not show any evidence for non- 
thermal emission which has been detected in the 
former system (Richards et al. 1999). 

The majority of optical papers in the literature 
have concentrated on spectroscopic observations 
and there is little in the way of imaging. What 
there is has revealed an extended nebula (Bang 
et al. 1992, Corradi et al. 1999) but no counter- 



Table 2: Pixel ranges used in the images shown 
in Figs. 1 and 2. Flux units for all filters 
are erg/s/cm 2 /A/arcsec 2 ; the F502N/F437N ra- 
tio map is dimensionless. 



Filter a Max. Flux 



F218W 


5.3 x 10- 


-14 


4.5 x 10" 


13 


F437N 


3.6 x 10" 


-15 


4.8 x 10~ 


14 


F469N 


7.2 x 10" 


-15 


1.2 x 10~ 


13 


F502N 


2.5 x 10~ 


-15 


5.6 x 10- 


14 


F547M 


1.5 x 10- 


-16 


2.7 x 10~ 


15 


F656N 


1.2 x 10- 


-15 


2.1 x 10- 


-14 


Ratio 


(min. 


) 


287 





part in either the infrared or ultraviolet has been 
observed. The resolution of ground-based optical 
telescopes has been insufficient for any observa- 
tions of the inner nebula at sub-arcsecond scales 
to have taken place. 

In this paper we present ffi>T/WFPC2 images 
of V1016 Cyg with resolution of <0.1 arcsec. We 
combine our results with VLA radio imaging and 
archival ultraviolet HST/STIS observations which 
enable diagnostic tests to be applied in order to 
determine various physical parameters. 

2. Observations 

2.1. HST/WFPC2 

Our HST observations took place on 1999 
November 4 as part of GO programme 8330 on 
symbiotic stars. Over the course of two orbits im- 
ages were made in six different filters - F218W, 
F437N, F469N, F502N, F547M, F656N - details 
of which can be found in Table 1 and in Biretta 
et al. (1996). These filters were chosen in order 
to probe the hot and cool continuum components 
of the symbiotic system (F218W and F547M re- 
spectively, see Eyres et al. 2001 for assumptions 
and justification), as well as various emission lines 
from which conditions within the nebula could be 
determined. 

The calibrated images were obtained from the 
HST archive. Since there were two observations 
in each filter it was possible to perform cosmic ray 
removal. A second pair of shorter exposures was 
also taken in each filter; this time the second image 
was 'dithered' in order to allow full recovery of the 
spatial resolution (Biretta et al. 1996). The iraf 
task DRIZZLE was used to recombine the images. 
The resultant pixel size of the dithered images was 
0.02275 arcsec, compared with 0.0455 arcsec for 
the undithcrcd images. We use the dithered im- 
ages to obtain the results below. 

Each image has been de-reddened using the 
standard interstellar extinction curve of Whittct 
(1992) and published values of Rudy et al. (1990), 
who determined a value of E(B-V)=0.46±0.007 
from the hydrogen lines. (See also Section 3.4.) 

2.2. HST/STIS 

We have also obtained two STIS spectra from 
the public archive (Proposal 8188, Cycle 8, P.I. 
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Fig. 1 — HST images of V1016 Cyg in the F218W ('&ot continuum', top left), F469N (He n A4686, middle 
left), F547M ('cool continuum', bottom left), F437N (A4364, top right) F502N (AA4959, 5007, middle right) 
filters. The sixth plot (bottom right) is the ratio of filters F502N:F437N. Greyscale pixel ranges for these 
images are listed in Table 2. The same axes have been used for each of the plots to aid comparison; the 
black line across each r>lot corresponds to 2000 AU. assuming a distance of 2kx>c to V1016 Cve. Contours 



Table 1: A list of the filters, their properties and the exposure times used in our HST/WFPC2 observations, 
f refers to those images which were obtained in the dithered mode in order to improve the spatial resolution. 



Filter 


Exposure 


A 


AA 


Peak A 


Scientific features 




time (s) 


(A) 


(A) 


(A) 


and wavelengths (A) 


F218W T 


40 


2136 


355.9 


2091 


Hot continuum 


F218W 


100 


2136 


355.9 


2091 


Hot continuum 


F437Nt 


40 


4369 


25.2 


4368 


[O in] A4364 


F437N 


510 


4369 


25.2 


4368 


[O in] A4364 


F469Nt 


20 


4695 


24.9 


4699 


He ii A4686 


F469N 


100 


4695 


24.9 


4699 


He ii A4686 


F502Nt 


20 


5012 


26.8 


5009 


[O in] AA4959,5007 


F502N 


100 


5012 


26.8 


5009 


[O in] AA4959,5007 


F547Mt 


20 


5454 


486.6 


5362 


Cool continuum 


F656Nt 


10 


6562 


22.0 


6561 


Ha A6563 




Fig. 2.— HST Ha image taken with the F656N filter (A6563, left) and the 23 GHz radio map (right). 
The greyscale pixel range for the Ha image is listed in Table 2; the approximate position of the slit is 
indicated (the position angle of which is coincidentally aligned roughly with that of the binary), as well as 
the extraction positions of the spectra shown in Fig. 4 (see Section 3.5). Contours on the radio map indicate 
radio flux density values of 3, 6, 12, 24, 48 and 96 xcr where a = 0.19 mJy. The black horizontal lines 
correspond to 2000 AU and 200 AU in the Ha image and radio map respectively, assuming a distance of 
2kpc to V1016 Cyg. 
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CD. Keyes). V1016 Cyg was observed on 2000 
March 24 using the gratings G140L and G230L; 
the wavelength coverage is 1145-1740 and 1570- 
3180 A with a spectral resolution of 0.60 and 1.58 
A/pixel respectively. The slit width and length 
were 0.5 and 52 arcsec respectively with a pixel 
scale of 0.025 arcsec/pixel on the CCD chip. 

2.3. VLA 

The VLA observations were made on 1999 
September 26 at 23 GHz. Whilst not strictly 
simultaneous with the WFPC2 observations we 
note that a ~ 40 day period is short, relative to 
the likely timescales of variability in the source, 
and so therefore it is possible to compare them 
directly. 

The images were reduced using standard flag- 
ging, calibration and imaging techniques within 
MPS. The primary calibrator was 3C286 (1331+305) 
which led to determination of the flux of the sec- 
ondary calibrator, 2007+404, as 1.57 ± 0.07 Jy. 
The complex gain solutions obtained were then 
applied to the target, V1016 Cyg for which a peak 
flux of 23.28 ± 0.17 mJy/beam was measured. 
The final image was naturally weighted in order 
to improve the somewhat high noise levels at this 
frequency 

3. Results and discussion 
3.1. Images 

The resulting dereddened images can be seen in 
Figures 1 and 2. Fig. 1 shows the HST images in 
the F218W, F469N (He n A4686), F547M, F437N 
([Om] A4364) and F502N ([Om] AA4959, 5007), as 
well as the ratio of the two [Om] images. Fig. 2 
shows the image taken with the Ha filter, F656N 
and the 23 GHz radio map. North is to the top of 
the image in each case and an AU scale is given 
for each image at an assumed distance of 2kpc (the 
black line represents 2000 AU in each HST image 
and 200 AU in the radio and map). 

The majority of the HST images show an ap- 
proximately circular (spherical?) nebula ranging 
in size from 0.5-2 arcseconds in diameter; there are 
very few additional features. However in the case 
of F502N (Fig. 1) the nebula is more asymmetric 
and, as we show in the next section, extended ap- 
proximately perpendicular to the derived binary 



axis. We must consider the possibility that this 
is some form of bipolar outflow. Closer inspection 
of the brighter regions of the image reveals more 
intricate features with various filaments and en- 
hanced 'blobs'. Comparison with the larger-scale 
images of Corradi et al. (1999) is difficult due to 
the different resolution scales - nonetheless, pos- 
sible extension in the same direction as their work 
is suggested. 

The 23 GHz radio map (Fig. 2) shows a source 
with peak intensity 23.3 ± 0.17 mJy/beam. Con- 
tours indicate flux density values of 3, 6, 12, 24, 
48 and 96 xa where a = 0.19 mJy Structure 
is observed on a smaller scale compared with the 
optical; we note that the spatial resolution of the 
VLA and dithered HST images are comparable. 
Again an extension approximately perpendicular 
to the binary position angle (see below) is ap- 
parent, as well as an additional extension rotated 
~ 90° clockwise from this. 

The ratio map and the annotations on the Ha 
image are explained in Sections 3.3 and 3.5. 

3.2. Stellar positions 

By making use of the dithering technique in or- 
der to recover full spatial resolution, we have been 
able to measure the position of the two stars of the 
symbiotic system. We use the argument presented 
in Eyres et al. (2001) to show that it is reasonable 
to assume that the hot component dominates the 
emission peak in the F218W filter and that the 
cool component dominates in the F547M filter. 

We have used the AIPS task imfit to fit a 
Gaussian to these two images - using of course 
the images in which the pixel size was 0.02275 arc- 
sec. We note that in each case the Gaussian pro- 
vided a good (xt ~ 1) fit to the radial profiles and 
was in agreement (within the errors) with com- 
parison point spread functions produced by the 
tinytim software. The peak position was deter- 
mined relative to 19 57 05 +39 49 36 - the shift 
in the F218W image was a = 0.09860 ± 0.00004, 
S = 0.9130 ± 0.00005; in the F547M image it was 
a = 0.10028 + 0.00004, 5 = 0.8789 + 0.00005. Thus 
the hot component and the cool component are 
separated by 42.4 mas (to within an accuracy of 
< 1 mas) at a position angle of 143.5° (to within 
an accuracy of < 10°). 

This is only the second time the separation of 
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the two stars in a symbiotic binary has been mea- 
sured directly (HM Sge being the first - Eyres et 
al. 2001) and such a measurement can place im- 
portant constraints on estimates of the distance 
to and orbital period of the system. We used Ke- 
pler's third law to calculate that for a projected 
separation of 42 mas, an assumed distance of 2 
kpc and an assumed total mass of 2M the min- 
imum orbital period is ~ 544 years. Since the 
system is non-eclipsing, thus ruling out a high in- 
clination angle and increasing the separation, the 
period will actually be higher than this; obviously 
using one of the higher published distance esti- 
mates increases the period further still. In order 
to decrease the orbital period to ~ 80 years (Schild 
& Schmid 1996) it is necessary to decrease the dis- 
tance to ~ 580 pc. We note that even a period of 
80 years is high compared with other symbiotic 
systems (Belczyhski et al. 2000) and suggest that 
either the distance to V1016 Cyg needs revising or 
the period of V1016 Cyg is unusually high or the 
published orbital periods for other symbiotic sys- 
tems may be too low. However, we also note that 
these published values tend to be for S-type sym- 
biotics which are expected to have shorter obital 
periods; D-type symbiotics are indeed expected to 
have periods on timescales of hundreds of years 
and so our result is consistent with this. A final 
possibility is that our HST observations took place 
around apastron of a highly eccentric orbit. 

At this point we must add a note of caution to 
these results. The values for angular separation 
and position angle of the two components in V1016 
Cyg and HM Sge are very similar (V1016 Cyg - 
42.4 mas, 143.5°; HM Sge - 40 mas, 130°), perhaps 
suggesting some, as yet unknown, instrumental ef- 
fect. We note that the -ffSTroll angle was similar 
for the two targets (252° vs. 259°). However mea- 
surement of the angular separation of the stellar 
components in CH Cyg (Eyres et al. 2002) has also 
been attempted and none was found, to within the 
errors of the method. This was consistent with the 
fact that the hot component was eclipsed by the 
giant at the time of the HST imaging, as shown 
by photometric monitoring, and adds credibility 
to the method. Furthermore, the position angle 
for the components of HM Sge is confirmed using 
a second method (Schmid et al. 2000). Perhaps 
we should not be too surprised by our results for 
V1016 Cyg as the observed similarities between 



them (e.g. position and separation of radio peaks; 
Watson et al. 2000, Eyres et al. 2001) would lead 
us to expect similar parameters. However further 
investigation of any possible instrumental cause is 
probably warranted. 

3.3. Nebular diagnostics 

The ratio of the [Oiii] lines in the nebula 
{Rom) depends on both electron temperature 
T e (K) and density n e (cm~ 3 ) (Osterbrock 1989) 
according to: 

R = JA4959 + JA5007 = 7.73ea:p[3.29 x lO 4 ^- 1 ] 
OIH ' JA4363 ""' 1 + 4.5 x 10- 4 (n e T e - a5 ) 

(1) 

Furthermore, the radio brightness temperature 
T b in Kelvins is given by 

T b = T e (l-e- T ") (2) 

where the optical depth r„ at frequency v 
(GHz) is 

t v w 8.24 x lO" 2 ^- 1 - 35 !/- 21 ^/ (3) 

where I (pc) is the path length through the neb- 
ula. Thus the three relationships can be used to 
plot the ratio Rom and the brightness tempera- 
ture Tb as functions of T e and n e . 

The ratio of the two [Oiii] maps has been ob- 
tained and is shown in Fig. 1. We note that 
the [Oiii] lines dominate the emission in the filter 
bands used for our observations (e.g. Andrillat, 
Ciatti & Swings 1982). Contours indicate Rom ~ 
2, 4, 6, 8 and increase outwards from the centre 
of the nebula. In the outer regions values increase 
further up to Rom ~ 15 and the occasional pixel 
has a value higher still, although these should be 
treated with caution. The higher value contours 
have been omitted from the plot for clarity. There- 
fore, by using the same assumptions and consider- 
ations as Eyres et al. (2001) we can use our ratio 
map to determine the conditions of the nebula at 
different path lengths. 

Curves for Rom =2, 4, 8, 16 have been plotted 
in T e :logn e space and can be seen in Fig. 3. The 
brightness temperature of the radio emission has 
also been determined for the second-highest (48cr, 
i.e. T b ~ 5424 K) and lowest (3er i.e. T b ~ 339 
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K) contours in the radio map. The 23 GHz radio 
image suggests that the angular size of the nebula 
is 0.1-0.5 arcsec, corresponding to 0.001-0.005 pc 
(or <~ 200-1030 AU assuming a distance of 2 kpc); 
unfortunately it is not possible to probe the more 
extended regions of the nebula without lower fre- 
quency and/or resolution radio images). Thus a 
further four curves corresponding to these values 
of Tf, and I have also been included in Fig. 3. 

Fig. 3 suggests that n e increases with increased 
brightness temperature and decreases with in- 
creased path length, as we would expect. The 
regions of lowest Rom are found right at the cen- 
tre of the source (see ratio map in Fig. 1) and it 
is clear that the regions of the nebula probed by 
our radio observations correspond to Rom < 8 
and electron temperatures in excess of 40-50000 
K - suggesting some shock heating may begin to 
become important here. The extended nebula is 
cooler, with Rom > 8; however, radio imaging 
at lower frequencies is required to probe these re- 
gions. We note that our results are not inconsis- 
tent with those of Schmid & Schild (1990) who de- 
termine a similar value of Rom but a temperature 
of ~ 30000 K. Their observations were at consid- 
erably lower spatial resolution and thus included 
the cooler regions of the nebula; our high reso- 
lution observations probe only the innermost and 
potentially significantly hotter parts of the nebula. 

3.4. Extinction 

We make use of the fact that the dereddened 
H/3 flux can be derived from the radio flux via the 
equation (Pottasch 1983): 



S v = 2.51 x 10 8 T e °- 53 ^- 01 F(H/3) Jy (4) 

Thus from our measured values of the Ha flux 
we can extrapolate the value for the measured H/3 
flux (Osterbrock 1989) and determine an estimate 
for the extinction across the nebula from the ratio 
of the radio to Ha images: 



E(B - V) 



1.46 



log 



S v 



F obs (Ha) 



- 7.15 (5) 



We note that, due to the poor absolute posi- 
tional accuracy of the HST compared with the 



VLA, it was necessary to assume that the peak 
flux position of both the radio map and Ha im- 
age were identical and to re-align the HST image 
accordingly. We also assumed that the radio emis- 
sion is purely thermal and optically thin - these 
are reasonable assumptions as outlined in Eyres et 
al. (2001). 

The resultant map showed extinction peaks ei- 
ther side of the radio peak as one might expect 
from a hot source surrounded by an extended 
dusty nebula. However the values of E(B— V) de- 
rived (ranging from 1.4±1.0 to 2.5±2.0) are high 
compared with previous estimates (e.g. Birriel et 
al. 2000), possibly due to self-absorption in Ha, 
and noisy due to the poor S /N in the radio image 
and the broad PSF of the Ha peak. 

3.5. STIS spectra 

The slit of the spectrograph lay at a position 
angle of 119.28° (from north measured through 
west) and had dimensions 0.5 x 52 arcsec; the ap- 
proximate position of the slit is indicated in Fig. 2. 
The zero position of the pointing was a = 19 57 
05.0, 6 =39 49 36.5, corresponding to pixels y — 
392 and 512 for the G140L and G230L gratings 
respectively. Thus we were able to determine that 
the slit lay roughly parallel to the derived binary 
axis but that it lay across the nebula, not actually 
over either of the two stars. 

Using the iraf package xId we extracted one- 
dimensional spectra at many positions along the 
slit; the intervals were determined by the signal- 
to-noise possible at each position, with a greater 
degree of averaging taking place towards the outer 
parts of the slit, particularly for the G230L grat- 
ing. Spectra centred on position y = 392 and 512 
(i.e. the slit centre) for the G140L and G230L 
gratings respectively) are shown in the top pan- 
els of Fig. 4 (averaged over 0.025 and 0.075 arc- 
sec along the slit respectively) ; the bottom panels 
show sample spectra 1.45 arcsec and 0.425 arc- 
sec from the slit centre respectively (averaged over 
0.12 and 0.25 arcsec along the slit respectively). 
Each of these three positions is indicated by dotted 
lines on the Ha image in Fig. 2. We then used the 
SPLOT task to measure the integrated flux under 
various emission lines for each extracted spectrum 
and compared these values across the nebula. In 
order to do this direct comparison it was necessary 
to take into account the averaging of rows during 
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Central 0.075" of G230L slit 



800 1500 2000 2500 



3000 3500 



Cut positioned at 1.45" (width 0.125") 
from centre of G140L slit 




C IV 
I 



. — ;1L 





1500 



Cut positioned at 0.425" (width 0.25"); 

from centre of G230L slit : 




Fig. 4. — Spectra from the centre of the G140L (top left) and G230L gratings (rows y = 392 and 512 respec- 
tively). These correspond to the brighter regions of the nebula which were probed during the observations. 
The two bottom panels show sample spectra from further out in the nebula - from rows y = 450 and 495 
for the G140L (bottom left) and G230L (bottom right) gratings respectively. These pixel values correspond 
to 1.45 arcsec (averaged over 0.125 arcsec) and 0.425 arcsec (averaged over 0.25 arcsec) from the slit centre 
respectively. The extraction positions are indicated by dashed lines on the Ha image in Fig. 2. 
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the extraction. Thus each flux was divided by the 
number of rows averaged, the pixel scale and the 
slit width, resulting in fluxes per square arcsecond. 

The resulting spectra from both gratings, al- 
though particularly G140L, showed that a num- 
ber of emission lines remained visible in the outer 
parts of the nebula - we show the integrated fluxes 
of lines Cn] A2326, Cm] A1909 and Civ A1549 in 
Fig. 5. As we would expect, the fluxes decreased 
from the centre of the slit outwards to the edge 
of the nebula. With a spatial scale of 0.025 arc- 
sec/pixel it is clear that the nebula maintains a 
significant flux level to a diameter of at least 15 
arcsec and possibly as much as 20-25 arcsec, al- 
though improved sensitivity is required to confirm 
this. This may seem surprising when compared 
with the images presented earlier in this paper 
which span about 2 arcsec at most. However we 
draw attention to work by Bang et al. (1992) who 
found evidence for a ~ 24 arcsec nebula from long- 
slit high-resolution spectroscopy of the Ha, He n 
A4686 and [O in] A5007 emission lines. The pres- 
ence of this extended emission was confirmed by 
Corradi et al. (1999) who obtained [N n] A6583 
spectra and images of V1016 Cyg; in this case the 
nebula appeared extended in the same direction 
as the [O ill] AA4959,5007 image above. This is 
the first time that the ultraviolet counterpart to 
the extended nebula has been discovered, although 
with only one STIS observation at a single slit 
position we are unable to investigate the detailed 
morphology of this extended ultraviolet emission. 

Finally we have plotted in Fig. 6 the Cn] A2326 
/ Cm] A1909 vs. Civ A1549 / Cm] A1909 ratios 
as for the diagnostic diagram of Allen et al. (1998, 
although we note that this model has been devel- 
oped for active galaxies). Points from the outer- 
most parts of the nebula have not been included 
due to the extremely high uncertainties in these 
regions. All other points lie in a region of the dia- 
gram which is suggestive of photoionisation, al- 
though the possibility of shocks is not actually 
ruled out. Perhaps surprisingly, there is no re- 
lationship between position in this diagram and 
position in the nebula, although it is possible that 
this will change with improved sensitivity data. 
These results are preliminary and more detailed 
modelling will take place, particularly if we ob- 
tain data with higher signal-to-noisc at both opti- 
cal and ultraviolet wavelengths. 



3.6. Comparison with HM Sge 

HM Sge and V1016 Cyg share a number of 
common properties which has led frequently to 
their simultaneous study (including e.g. Corradi 
et al. 1999, Schmid & Schild 1990). Both sys- 
tems are D-type symbiotics with a Mira as the 
cool component; furthermore they are both sym- 
biotic novae which have experienced an outburst 
within the last forty years. They are both rela- 
tively bright radio sources and the published es- 
timates of their orbital periods are relatively long 
compared with other symbiotic systems. Compar- 
ison of our results above with Eyres et al. (2001) 
suggests that coincidentally both central binary 
systems had similar separations and position an- 
gles at the time of our observations. 

Unfortunately such comparisons arc insuffi- 
cient if we wish to study the characterstics of 
the sources. Despite such similar parameters HM 
Sge is considerably brighter and appears to have 
much more structure within the nebula. Fur- 
thermore there is evidence for non-thermal radio 
emission from extended regions whereas there is no 
such emission in V1016 Cyg observed to date. In 
addition, our preliminary ultraviolet diagnostics 
suggest that photoionisation is more significant 
than shocked emission in the V1016 Cyg nebula. 
The ~ 10-year gap in time between the oubursts 
of these two sources may explain some of these 
differences. Current distance estimates to V1016 
Cyg (2-10 kpc) suggest that it is considerably 
further away than HM Sge (~ 1250 pc, Richards 
et al. 1999) and this may explain some of the 
differences; however if, as we suggest above, the 
distance is actually less than those published then 
it would be interesting to determine more about 
how and why these two sources differ, rather than 
focussing soley on their similarities. 

4. Conclusions 

We have presented the first sub-arcsecond opti- 
cal images of the symbiotic star V1016 Cyg and 
have shown that although the source was ap- 
proximately circular in most images, the [O III] 
AA4959,5007 emission appears to trace a bipolar 
outflow with axis aligned approximately NE-SW. 
These images have allowed us to determine the 
positions of the two stellar components of the sys- 
tem and thus a separation of 84 AU (assuming a 
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distance of 2 kpc) and position angle of 143.5°; 
the subsequent orbital period calculated for such 
a projected separation is extremely long (~ 544 
years) and may require revision of the distance es- 
timate and/or a highly eccentric orbit. We have 
added a note of caution here, however, as we 
need to conduct further work to totally eliminate 
any possible instrumental effects. Furthermore we 
have presented a 23 GHz radio image which we 
have used in calculation of firstly the electron den- 
sity and temperatures for different path lengths 
through the nebula and secondly the extinction at 
different locations within the nebula; the bright- 
ness temperature near the centre of the radio im- 
age was 5424 ± 113 K. Finally we have determined 
accurately the position of the slit in archival STIS 
spectra thus enabling us to show that there is ex- 
tended emission in the ultraviolet ina^ 15 — 25 
arcsec (diameter) nebula, as observed in optical 
data previously and which appears to be pho- 
toionised from preliminary analysis. Further STIS 
spectrosocopy, particularly if used to 'map' the 
nebula, is strongly encouraged in this and other 
extended nebulae of symbiotic sources. 
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Fig. 3. — Plots of logn e against T e in terms of 
Ron i (solid lines) and Tb, I (dotted/dashed lines). 
The different values of Rom, Tb and I are indi- 
cated next to the curves and in the key (see text 
for details). 
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Fig. 5.— Integrated fluxes for the Cn] A2326, Cm] 
A 1909 and Civ A 1549 emission lines plotted as a 
function of radial distance from the slit centre. 
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Fig. 6.— Ratios Cn] A2326 / Cm] A1909 and Civ 
A1549 / Cm] A1909 plotted as for diagnostic dia- 
gram of Allen ct al. (1998). With the exception 
of points to the edge of the nebula (for which the 
error-bars are extremely large) all points appear 
to lie in the region of photoionisation from their 
models, although without totally ruling out the 
possibility of shocks. 
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